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HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China’s space 
station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy 
resolution (~ 1% for electrons and photons and 20% for nuclei) and a large geometry factor (> 3m^ sr for 
electrons and diffuse photons and > 2 m^ sr for nuclei). In this work we discuss the capability of HERD to 
detect monochromatic y-ray lines, based on simulations of the detector performance. It is shown that HERD 
will be one of the most sensitive instruments for monochromatic y-ray searches at energies between ~ 10 to 
a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large 
effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric 
dark matter model can be probed with HERD. 

PACS numbers: 95.35.+d,95.85.Pw 


I. INTRODUCTION 


Cosmological observations have well established that dark 
matter (DM) constitutes ~ 25% of the Universe’s enemy con¬ 
tent and dominates the ordinary, baryonic matter |[ll|2l- The 
search for the DM particle becomes one of the most important 
tasks in the modern physics. The high energy monochromatic 
y-ray emission would be a “smoking gun” signature of parti¬ 
cle DM dH] . With the remarkable improvement of the sensitiv¬ 
ity and energy resolution of the y-ray detection by space and 
ground-based y-ray instruments, great efforts have been un¬ 
dertaken to search for monochromatic y-rays or sharp spectral 
features in the past few years l&fil. There is no compelling 
evidence for the existence of line-like y-rays yet. 

Increasing the energy resolution is very crucial for the 
monochromatic y-ray detection. The energy resolution of the 
current y-ray detectors, such as the Fermi Large Area Tele¬ 
scope (Fermi-LAT) in space and the Imaging Atmospheric 
Cherenkov Telescope (lACT) arrays on the ground, is of the 
order of 10% - 15% for (9(100) GeV photons. Such a reso¬ 
lution is not enough to firmly identify a y-ray line, when the 
photon statistics is not very high |la,0- The Alpha Magnetic 
Spectrometer (AMS-02) onboard the International Space Sta¬ 
tion has an energy resolution of ~ 2% at (9(100) GeV, but the 
effective area of AMS-02 is too small to search for the weak 
signals 01. The next generation of space-borne high energy 
cosmic ray (CR) and y-ray detectors, including the CALori¬ 
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metric Electron Telescope (CALET)' lfH . the DArk Matter 
Particle Explorer (DAMPE)2 01 and GAMMA-TOO^ 01 are 
designed to perform very high energy resolution (~ 1% - 2%) 
detection of photons with large effective areas. On the other 
hand, at TeV energies the ground-based Cherenkov Telescope 
Array (CTA, 01 ) will improve the capability of line-like y- 
ray searches. See Refs. 0H3 for the expected performance 
of these future experiments. 

The High Energy cosmic-Radiation Detection (HERD) fa¬ 
cility onboard China’s space station has been proposed re¬ 
cently ll2l . HERD is basically a five-side active calorimeter 
designed to perform high energy resolution and high statistics 
measurements of the CR nuclei, electrons and positrons, and 
y-ray photons in space. The scientific objectives of HERD in¬ 
clude the high sensitivity search for particle DM, direct mea¬ 
surements of the CR nuclei spectra and composition up to 
knee energies, and high energy y-ray sky surveys. Based on 
the detailed simulations of the HERD detector 11^ . we inves¬ 
tigate the expected performance of HERD on the monochro¬ 
matic y-ray line detection and the potential on the constraints 
of DM model parameters in this work. 

This paper is organized as follows. In Sec. II we briefly 
introduce the design and performance of the HERD detector. 
The sensitivity of line searches is presented in Sec. III. Tak¬ 
ing the Minimal Supersymmetric Standard Model (MSSM) 
as an example, we show the capability of HERD to explore 
the corresponding DM parameter space in Sec. IV. Einally we 
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conclude in Sec. V. 


II. HERD DESIGN AND PEREORMANCE 

HERD is composed of a 3-D cubic calorimeter (CALO) 
surrounded by microstrip silicon trackers (STKs) from five 
sides, while the bottom is left for mechanical support. Fig. [T] 
shows a schematic plot of the basic design and the major func¬ 
tions of each part of the detector. The CALO is made of 
21 x21x21 cubicLYSO crystals with each cell 3x3x3 cm^ 
which is coupled with the wavelength shifter fibers and read 
out by ICCD. This design corresponds to about 55 radiation 
lengths and 3 nuclear interaction lengths, respectively. The 
CALO detects the total energy deposited by the electromag¬ 
netic (EM) shower from CR leptons/photons and the hadronic 
shower from CR nuclei, and provides lepton/hadron separa¬ 
tion through the differences of the shower shape. The five 
sides of identical STKs, each consists of seven layers laid per¬ 
pendicularly one by one, are sandwiched with tungsten con¬ 
verters (2 radiation lengths), to ensure the maximum field of 
view (FOV) and provide charge identification, trajectory mea¬ 
surements, back scattering rejection, as well as some early 
shower information of photons and electrons. The detector is 
surrounded by plastic scintillators from five sides, which are 
used to reject the low energy charged particles to maximize 
the efficiency for photons and high energy CRs including elec¬ 
trons. 



FIG. 1: Schematic design plot of HERD. 


Extensive simulations have been carried out with GEANT4 
and FLUKA to get the scientific performance of 
HERD Q. We generate events, including photons, electrons 
and protons, at energy grids of 0.5, 1, 5, 10, 50, 100, 200, 300, 
500, 700, and 1000 GeV, with an isotropic event generator. 
The shower development in the detector of each event is simu¬ 
lated using the Quark-Gluon String Precompound (QGSP) in¬ 
teraction model for photons/electrons and DPMJET-III model 
for protons. The size, orientation, and shape of the shower 
is used to reconstruct the energy, direction, and type of the 


incident particle. To ensure the best determination of the en¬ 
ergies for photons and electrons, we select those events whose 
shower maximum is fully sampled in CALO, which results in 
a decreasing fiducial area with increasing energy. The lepton- 
hadron discrimination is primarily based on the shower shape 
in CALO. The charge measurement by STKs can be used 
to identify neutral particles from charged ones. A machine 
learning, multivariate analysis with boosted decision trees is 
adopted to identify the particle type. The background (elec¬ 
trons and protons) rejection efficiency is obtained when keep¬ 
ing 90% signal (photons). We assume infinite shield of the 
bottom side of the detector from the space station. The shield 
of the other five sides is negligible. This assumption will ef¬ 
fectively exclude the earth limb’s photons due to the fact that 
the space station will point away from the earth surface for 
almost all of the time. The secondary events due to inter¬ 
actions between incident CRs and the space station, which 
would travel upwards in the detector, will also be excluded. 

The ideal energy resolution of the EM shower is about 0.1% 
at 200 GeV. Considering the stochastic fluctuation of the num¬ 
ber of photoelectrons and the error of energy calibration, the 
energy resolution for the EM component can be ~ 1%. Due to 
the large nuclear interaction length of the CALO, the energy 
resolution of the hadronic component is around 20%, almost 
constant from hundreds of GeV up to PeV. The CALO can be 
also used for the lepton/hadron separation due to the fact that 
the EM and hadronic showers differ in their spatial and energy 
distributions in the high granulated crystals. The hadron effi¬ 
ciency is about 5 x 10~® when keeping 90% of the EM events. 
The key performance of HERD, in comparison with previous 
and other proposed missions, is its extremely large effective 
geometry factor for all types of high energy particles, primar¬ 
ily due to its thick 3-D CALO and five-sided STKs. With 
homogeneous design for detecting particles from every un¬ 
blocked direction, the effective geometry factor is > 3 m^sr 
for electron and diffuse y-rays, and > 2 m^sr for CR nuclei. 

The energy resolution of CALO can be parameterized as the 
quadratic addition of the statistical fluctuation term (oc E^C^) 
and the electronic noise term (oc £ *) 

The minimum energy resolution is taken to be 1%, due to 
the systematic effect of the calibration and the non-uniformity. 
The effective exposure depends on the orbit of the space sta¬ 
tion. Employing the orbit of Tiangong-1 in 2012-2013, we 
simulate the all-sky exposure map of HERD, taking the shield 
of the space station into account. The Galactic center region 
is of particular interest for the DM searches. An average ex¬ 
posure of ~ 0.65 yr is achieved for this direction and one 
year of operation. The photon detection efficiency is the prod¬ 
uct of the pair conversion efficiency (~ 80%), the recon¬ 
struction efficiency (~ 90%), the efficiency for electron rejec¬ 
tion (90% for 10“^ residual electrons), and the efficiency for 
hadron rejection (90% for 10 * residual hadrons), which is in 
total ~ 60%. The geometry area of HERD/CALO is 63 x 63 
cm^ on each side. If the energy of the incident y-ray photon is 
high enough, there will be leakage of the shower for the events 
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hitting the edge of the detector. The fiducial area for good 
events can be approximated as (3660 - 717 x log(£’/GeV)) 
cm^, for the energies from ~ 10 GeV to ~TeV. Considering 
the correction of the shower leakage, the effective area can be 
larger. Here we adopt this result as a somehow conservative 
estimate. 


III. SENSITIVITY OF THE LINE SEARCH 
A. Backgrounds 

The backgrounds for monochromatic y-ray searches in¬ 
clude misidentified charged particles (nuclei and leptons) 
and the continuous y-ray events. The misclassified CRs are 
isotropic, due to the loss of the direction information during 
the diffusive propagation. For the proton flux, we adopt a fit¬ 
ting formula based on the recent AMS-02 data 12^ 

(j)p{E,) = 7.8 X 10-\£,/GeV)-0-® 

X [l.O -H (£i/4.9GeV)‘'^’]^‘ GeV^'cm^^s^^sr"'. 

( 2 ) 

The nucleon flux from Helium, is lower by a factor of ~ 3 - 5 
in the energy range of 10 - 1000 GeV/nucleon ll^ . Heav¬ 
ier nuclei play an even less important role compared to pro¬ 
tons. Since the hadron rejection power of HERD is very high, 
the nucleon background is always subdominant in the current 
study. The total spectra observed by AMS-02 and 

HESS m can be fitted by 

4>e4Ek) = 1.85 X 10^^(£*/GeV)^°-’‘ 

X [l.O + (Ek/3.5 

X [l -H(£i:/1300GeV)^]^°'^GeV^'cm^V'sr^'. 

(3) 

The efficiencies of protons as well as and e , rjp - 10 * and 
rje = 10“^, are multiplied with the CR fluxes in order to get 
the backgrounds of the misclassified CRs. With such a high 
rejection power of the CR hadrons, the hadronic background 
can be safely neglected for HERD. This is different from the 
lACTs, for which the hadronic background is comparable to 
other backgrounds S. 

The y-ray sky is composed of diffuse emission and point or 
extended sources. Since the bright point-like sources can be 
effectively removed in the data analysis, and the residual weak 
sources are expected to contribute < 10% to the diffuse events 
|@], we disregard the point sources in this analysis. Only the 
Galactic and extragalatic diffuse y-rays are taken into account 
as the continuous background. 

The Galactic diffuse y-ray emission can be modeled with 
three components that can be distinguished according to their 
radiation mechanisms; (a) the decay of neutral pions produced 
by inelastic collisions of CR protons with the interstellar gas, 
(b) inverse Compton scattering of the interstellar soft photons 
by CR electrons and positrons, and (c) bremsstrahlung radia¬ 
tion produced by the scattering of CR electrons and positrons 


with the interstellar gas. There are also a few identified large 
scale structures such as the Eermi bubbles ll^ 1^ . Loop I 
jm, and the Magellanic stream. The majority of the Fermi 
diffuse y-ray emission can be relatively well modelled with 
the CR interactions, the large scale structures, the residual 
CR events and point sources, except for some excesses in 
the Galactic plane ll^ . Based on the physical modeling of 
the Galactic diffuse y-ray emission and a likelihood fitting of 
the all-sky data, the Fermi Collaboration built diffuse y-ray 
templates for the point source analyses. Those templates are 
suitable to simulate the y-ray background in order to forecast 
the sensitivity of line searches with future experiments. Since 
the Fermi data are limited to a few hundred GeV, we need 
to extrapolate the templates to higher energies. We focus on 
the inner Galaxy region, where the uncertainties of the spec¬ 
tra of the Fermi bubbles will affect the extrapolations 11^ . 
Thus we adopt two different templates, to take such uncer¬ 
tainties into account. The first template is the one used for the 
“p6vll” data analysis^ without the Fermi bubbles. The sec¬ 
ond one is built with the Fermi bubbles and is used for the Pass 
7 Reprocessed Data^. These templates are transformed into 
HEALPix® projections with Aside = 64. Then the spectrum of 
each pixel is extrapolated to 2 TeV. Due to the suppression of 
the Klein-Nishina cross section, the simple extrapolation will 
tend to over-estimate the inverse Compton scattering emission 
at high energies, which makes the line limits more conserva¬ 
tive. 

For the extragalactic diffuse emission, we adopt the fitting 
formula of the latest Fermi data, which extends up to 820 GeV 

M 

~ ‘ (0.1 GeV) 

The parameters depend on the assumptions of the foreground 
model i^ . In this work we adopt /qj = 0.95 x 10 
GeV^'cm^^s^'sr"', y = 2.32 and Ecut = 279 GeV, corre¬ 
sponding to model A in Ref. j^ . Different choices of the 
parameters will not affect our results significantly because the 
extragalactic background is only a sub-dominant contribution 
to the total background. 

Fig. 12 shows the contributions of different components 
to the background. Since the Galactic diffuse emission is 
anisotropic, the spectra are given for different directions in 
the sky, concretely for the Galactic Center, as well as for low 
and high latitude regions. We can see that at low latitudes, the 
contribution from the Galactic diffuse component dominates 
over the other components. When approaching the high lat¬ 
itude regions, the misclassified and become more and 
more important, especially at low energies. The fluxes of mis¬ 
classified protons and the extragalactic y-ray background are 
always lower than the and backgrounds. 


^ http://fermi.gsfc.nasa.gOv/ssc/data/analysis/software/aux/glLiem_v02_P6_V 11 
_DIFFUSE.fit 

^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/aux/gll_iem_v05_revl.fit 
® http://healpix.jpl.nasa.gov 
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FIG. 2: Comparison of the different backgrounds; the Galactic and isotropic diffuse confinuous y-ray backgrounds, CR electrons, and protons 
(multiplied by 10^). The left panel is for “p6vll” template, and the right panel is for “p7v6” template. The Galactic diffuse y-ray emission is 
anisotropic, and we show the fluxes in three directions with (/, b) = (0,0), (0,20°), and (0,40°), respectively. 


B. DM profile and region of interest 


The density profile of the DM distribution in the Milky Way 
halo comprises large uncertainties in the innermost region. As 
generally adopted in the literature, we consider several typical 
DM density profiles to take such uncertainties into account. 
The first one is the commonly adopted cuspy Navarro-Frenk- 
White (NFW) profile il, 


p(r) = 


Ps 


(5) 


with r, = 20 kpc. The second one is the Einasto profile 
with an asymptotic flat slope in the center 


p{r) = Ps exp{-(2/a)[(r/rX - 1]), (6) 


where = 20 kpc and a - 0.17, which is favored by more 
recent simulations BH. The third one is the cored isothermal 
profile ll4^ 


p{r) = 


Ps 

1 -I- (r/r^)2 


(7) 


with G = 5 kpc. All the profiles are normalized to local den¬ 
sity of po = 0.4 GeV cm“^. 

Optimized regions of interest (ROI) are very important in 
the search for the weak DM signal. Since the expected signal 
depends on the DM density profile, we adopt different ROIs 
for different assumed profiles. Following Ref. |@], we choose 
a circular region with radius of 16° (R16) for the Einasto pro¬ 
file, 41° (R41) for the NEW profile, and 90° (R90) for the 
isothermal profile, in the case of DM annihilation. We also 
discuss the decaying DM scenario, for which the optimized 
search region (R180) corresponds to the whole sky for all 
these profiles. In fact, the integrals of the DM density over 
the Milky Way are very similar among these profiles (see Ta¬ 
ble |I|i. Eor definiteness, we adopt the NEW profile for decay¬ 
ing DM. In all these ROIs the Galactic plane with |/| > 6°, 


\b\ < 5° is removed. Eig.|3]shows the ROIs adopted to search 
for monochromatic y-rays in this work. Table |T] summarizes 
the optimized ROI for each DM halo profile and the corre¬ 
sponding 7-factors (integral of p^ or p over the line-of-sight 
and ROI). The results differ slightly from that of Ref. |@1, due 
to the mask of point sources in Eermi’s paper. 


TABLE I: Summary of the optimized ROIs and 7-factors. 


Profile 

ROI 

J ann 

ROI 

•^dec 



(10^^ GeV^srcm-^) 


(10^^ GeVsrcm-2) 

Einasto 

RI6 

9.39 

RI80 

2.55 

NFW 

R41 

9.17 

RI80 

2.52 

Isothermal R90 

6.95 

RI80 

2.56 


Mollweide view 



FIG. 3: Illustration of the different ROIs for different assumptions of 
the DM density distribution. From inside to outside, the regions are 
R16, R41, R90, and R180, respectively. 
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C. Mock Data and statistical treatment 


Given the expected y-ray fluxes including the CR back¬ 
grounds and the performance of the instrument, we can gener¬ 
ate mock data for HERD. We set 400 energy bins per decade 
in the energy range 5 GeV - 2 TeV, and calculate the expected 
number of counts in each energy bin i with width AE, is 

= At dE J dE'R(E,E')A,ff(E')(l>i:i(E'), (8) 

where At is the exposure time, R(E, E') is the energy response 
function of the instrument, A^s(E') is the effective area, and 
O™ is the total y-ray and CR background in the ROI. In this 
work we adopt 5 years of survey time, which corresponds to 
an average ~ 3.25 year effective exposure. The energy re¬ 
sponse function is assumed to be Gaussian, and its width is 
given in Eq. ([TJ- Assuming background only, we generate the 
mock observational counts in each energy bin, by gen¬ 
erating random numbers that are drawn from a Poisson dis¬ 
tribution with expectation Here we neglect the effect of 
the point spread function (PSF), which is expected to be small 
since we integrate the y-ray fluxes in large enough sky regions 
compared to the resolution angle. 

Then we calculate the expected number of counts of the the¬ 
oretical model with the DM contribution. The monochromatic 
y-ray flux from DM annihilations or decays reads 


<<rr>xJanni ^ R(jfi £), annihilation 

Smii^ (9) 

^^xR(m^,E), decay 

in which (crv) or t is the annihilation cross section or decay 
lifetime of the DM particle, is the DM mass, and Ny is the 
multiplicity of one annihilation or decay. On the other hand, 
the theoretical background is parameterized as a single power- 
law function = CE '>', which is good enough to approx¬ 
imate the background spectrum in a narrow energy window 
(see below the definition of the energy window). Our model 
for the y-ray flux contains thus three parameters, the two back¬ 
ground parameters C and y, as well as the flux from DM anni¬ 
hilation or decay (proportional to (crv) or 1/t). The expected 
number of counts are calculated by substituting the back¬ 
ground with our three parameter flux model, -H in 
Eq. (|8j. The likelihood function is given by 




Xocfl 


(«ixp)"* exp(-«',,p) 


obs 


( 10 ) 


We adopt the profile likelihood analysis jdffll and implement 
the sliding window technique i^, in order to derive the upper 
limits of the DM flux. We follow Ref. to define the en¬ 
ergy windows. For a specified central energy E, which is the 
energy of the y-ray line to be searched, the energy window is 
defined to be {E / y/f, {E Ve] ■ The parameter e is determined 
to ensure that the astrophysical background is well described 
by a power law. Specifically, we employ fittings to the 300 
mock data sets with a power-law model, for each ROIs and for 


different window sizes (from e - 1.2 to e = 8). The resulting 
distribution of the 300 values for each ROIs is then tested 
against a x^ distribution with the number of degrees of free¬ 
dom given by the number of energy bins in the energy window 
minus 2 (number of fitting parameters). An energy window 
is rejected if the corresponding p-value is smaller than 0.01. 
And our chosen energy window corresponds to the largest 
window that fulfills that criterion for all ROIs under consid¬ 
eration. This procedure gives e ~ 1.6 at 20 GeV and ~ 4 at 
1 TeV. We have tested that adopting the parameter e twice as 
large as the above derived values, the limits improve by only 
7% above 40 GeV. A possible optimization is to choose differ¬ 
ent e for different ROIs. However, the results will not change 
significantly. 

Maximizing the logarithm of Eq. ( fTOl i with respect to the 
two background free parameters (C, y) and one DM parame¬ 
ter best-fit log-likelihood log Xbf ■ The one¬ 

sided 95% confidence level (CL) upper limit on the DM flux is 
found by increasing the signal starting from the best fit value 
until -21ogT:((I)™^) = -21ogi:bf + 2.71 JH. 


D. Results 

We perform a scan of the photon energies from 10 GeV to 
1 TeV, and calculate the profile likelihood with respect to the 
monochromatic line flux for each photon energy. The 95% CL 
upper limits of the line fluxes for different ROIs are shown in 
Fig.g] Here we adopt the “p7v6” Galactic background model. 
The results are the geometric mean of the limits obtained from 
300 mock data realizations. It is shown that at low energies the 
differences among various ROIs are larger than at high ener¬ 
gies. This reflects the effect of the relative weights between 
the electron background and the Galactic diffuse emission. At 
low energies the electron background plays an important role 
in the total background. Thus the backgrounds in different 
ROIs scale with the solid angles of the sky regions, and hence 
differ significantly. At high energies, however, the anisotropic 
Galactic diffuse background, which is dominated by low lat¬ 
itude emissions, becomes more and more important, and the 
differences among these ROIs become smaller. 

The upper limits of the line fluxes can be easily translated 
into the upper limits of the DM annihilation cross section or 
decay lifetime. The left panel of Fig. 0 shows the 95% up¬ 
per limits on the cross section of DM annihilation into a pair 
of photons (crv)yy. We find very stringent limits, reaching 
cross sections of 6 ■ 10“^*^ ~5 ■ 10“^^cm^s“* at 10 GeV, and 
4 ■ 10“^^ ~ 9 ■ 10“^^ cm^s^' at TeV. For the two Galactic back¬ 
ground templates we find very similar results in all three re¬ 
gions. The largest difference amounts to about 30%, due to 
the contribution of the Fermi bubbles to the background (see 
Fig. |2]i. We expect that below ~ 100 GeV the “p7v6” back¬ 
ground should describe the actual background better, while 
above ~ 100 GeV where the energy spectra of the Fermi bub¬ 
bles show a cutoff S, the extrapolation of the “p7v6” tem¬ 
plate will over-estimate the background. The limits for R16 
and R90 differ by a factor of about 5 at 10 GeV and 2 at 1 TeV, 
even though the 7-factors of these regions are very similar. 
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FIG. 4: 95% upper limits on the monochromatic photon fluxes for 
different ROIs, where the Galactic background template “p7v6” is 
adopted. 


Such differences come from the different background levels 
of these ROIs. In analogy to the flux limits, larger differences 
at low energies arise due to the electron background. 

In the right panel of Fig. |5] we compare the HERD limits 
with the results from other current or upcoming y-ray facili¬ 
ties. All the limits shown are for the Einasto profile. For the 
HERD limits we adopt the “p7v6” background model. The 
Fermi results are adopted from Ref. §], with the same ROI 
(R16) and the newest analysis of the Pass 8 data. The predic¬ 
tion of GAMMA-400 is from a 20° region around the Galactic 
center, excluding the Galactic disc (|/| > 5° and \b\ < 5°), for 
5 years of full sky surve vKdIl . It is shown that HERD can 
improve the Fermi limits HI from 5.8 years of observation by 
up to a factor of a few. The limits expected from GAMMA- 
400 and HERD are comparable. However, the effective area 
of GAMMA-400 used in Re f. as well as the y-ray effi¬ 
ciency, seems to be too ideal MH . The lACTs are expected to 
be more effective to probe the y-ray lines at high energies. We 
also show the limits from 112 h of Galactic center observation 
with the HESS telescopes id in a circular 1 ° region around 
the GC, excluding \b\ < 0.3°, as well as the expected limits 
for CTA 12^ . The CTA limits are derived for the same region 
and observation time as that of HESS and are properly 
rescaled to y-ray lines. We can see that below ~ 300 GeV 
HERD is more sensitive than CTA, whereas at higher ener¬ 
gies CTA is expected to be more sensitive. 

The constraints on the lifetime of DM decay into yv are 
shown in Fig.| 6 l Here we adopt the NFW profile and the R180 
region. For the sake of comparison with Fermi, we show the 
limits on Tyy, for x ^ 7"^ channel. A classical example of 
this kind of DM is the gravitino in the supersymmetric model 
id^ll . Similarly, we find that the HERD limits can improve the 
current Fermi-LAT constraints §] by a factor of ~ 2 - 3. 


IV. CONSTRAINTS ON MSSM PARAMETERS 

In this section we show the HERD potential on the con¬ 
straints of the specific MSSM DM model parameter space. 
Since we are only interested in the y-ray line features, we fo¬ 
cus on the lightest neutralino in the MSSM, the DM candidate, 
that can monochromatically annihilate to yy and yZ via loop 
processes. Because of the loop suppression, the cross section 
is usually rather small compared to the total annihilation cross 
section. On the other hand, we want to demonstrate the power 
of HERD in the higher cross section region. Therefore we 
ignore the sfermion-DM coannihilation region that exhibits 
low annihilation cross sections. The mass splitting between 
a sfermion and DM is small enough to release the Boltzmann 
suppression so that the correct relic density can be reached by 
accounting for neutralino-sfermion coannihilation processes. 

In this analysis we consider the DM sector of the MSSM, 
where we take into account the bino M\, wino M 2 , and 
higgsino fi mass parameters. The p-term is assumed to 
be positive. Furthermore, the LHC multijet plus miss¬ 
ing energy search can put a very stringent mass limit 
on the gluino mass M 3 and the squark masses. To avoid 
this limit, we take the gluino mass to be the same as 
the universal sfermion masses my that are heavier than 
max[Mi, Ma,//, niA, 800 GeV]. The pseudo-scalar Higgs mass 
niA and the value of tan p control the higgsino mixing and the 
A-resonance region which can have higher annihilation cross 
sections into monochromatic y-rays. We unify all the trilin- 
ear couplings to be Aq. Finally, we show the ranges of the 
MSSM7 parameters used in our scan: 

3 < tanjS < 62, 

10 < (Mi,M 2 ,p)/GeV < 4000, 

200 < Ota/G eV < 8000, 

max[Mi,M2,jU,OT/i,800.0GeV] < mj < 8000GeV, 

M3 = mj, 

-5<Ao/TeV<5. (11) 

We assume a non-thermal relic scenario in which the neu¬ 
tralino can be reproduced by late decays from other parti¬ 
cles. The advantage of this scenario is to allow for a smaller 
relic density during the thermal freeze-out stage, and hence 
a larger interaction cross section of DM. The late time pro¬ 
duction of DM can be responsible to the measured relic den¬ 
sity today. The upper limit of relic density is taken from the 
PLANCK measurement H. Besides the relic density, we 
also consider some other common constraints. We only al¬ 
low chargino masses greater than 103.5 GeV The Higgs 
mass is constrained using the most recent combined analy¬ 
sis from ATLAS and CMS ll^ taking into account 2 GeV 
theoretical uncertainties in the Higgs likelihood function. In 
the neutralino mass region close to the ZJh mass, we con¬ 
sider the Z and Higgs invisible decays Moreover, one 
can the use the measurement of ^ b ^ sy 

to further constrain the Higgs sector. Here we use the up¬ 
dated data BR(b ^ sy) x 10^ = 3.43 + 0.22 + 0.21 jUl and 
BR(Bs ^ juV) X 10'^ = 2.9 + 0.7 HI. In the B, ^ 
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FIG. 5: Left panel: the mean limits obtained from 300 realizations of mock data, for the three ROIs (DM prohles): R16 (Einasto), R41 (NFW) 
and R90 (isothermal). The solid (dashed) lines correspond to “p7v6” (“p6vll”) Galactic background. Right panel: comparison of HERD 
results with limits from Fermi (solid is the observational limits and dashed is the expected mean) and HESS in, and that predicted for 
GAMMA-400 and CTA |2^ . See the text for details. 



FIG. 6: 95% upper limits of the decay lifetime Tyy of the DM parti¬ 
cle. The observed and expected mean limits from 5.8 years of Fermi 
observations are also shown for comparison ItIi . 


likelihood function we also include 10% theoretical uncer¬ 
tainties. Regarding DM direct detection which can put a 
stringent limit on the DM parameter space, we use the spin- 
independent cross section limit from LUX ifs^ . the DM- 
proton spin-dependent cross section limit from PICO-2L if^ 
and the DM-neutron spin-dependent cross section limit from 
XENON 100 il. 

The scan of the parameter space is performed using 
the package MultiNest ifs^ . which is optimized for the 
Bayesian sampling. The scans are driven by the likelihood 
function with all the above constraints. We use 15000 live 
points for the sampling. The evidence tolerance factor is 10“^, 
and the sampling efficiency is 0.8. To obtain a good coverage 
of the parameter space, we combine six separate scans, three 


of which have log priors of the mass parameters and the other 
three have flat priors. The SUSY mass spectrum is computed 
by using SOFTSUSY if^ and then passed to DarkSUSY ifs^ for 
DM observables, Superlso 11^ for Bs —> and b ^ sy 

computation, and SUSY-HIT ll^n for Higgs decay width. In 
addition, we have checked all our collected points under the 
package HiggsBounds S to ensure that we do not violate 
any existing Higgs constraints. 

It is easier to identify the neutralino features by consider¬ 
ing bino-like, wino-like, higgsino-like and mixed neutralinos. 
The neutralino is decomposed into bino, wino, and higgsi- 
nos as 

= XbinoB + ZwinoW + Zh„Hu + ■ 

The fraction of each component, g, where i denotes bino, 
wino, or higgsino, is dehned as gbino = gwino = 
and ghiggsino = Zjj + Zjj^. We, therefore, identify the neu- 
tralino composition as bino-, wino- or higgsino-like when the 
corresponding fraction g, > 0.95. The other linear combina¬ 
tions are denoted mixed neutralinos. 

In Fig. 121 we show the 2cr allowed scattering points for all 
the constraints described above, where for the relic density 
the 2cr upper limit is applied. The red, green, and purple lines 
are the HERD sensitivity for the R16, R41, and R90 regions 
(see Fig. [3]). Note that the Z-boson and SM Higgs boson res¬ 
onance do not appear in both plots since they have v^ low 
yy and yZ annihilation cross sections in the MSSM In 
the neutralino mass region between the Higgs resonance re¬ 
gion {m^ ^ mhl2) and the chargino-neutralino coannihilation 
region ^ m^± > 100.0 GeV, the DM relic density is over 
produced due to the absence of reduction mechanisms for the 
annihilation cross section at freeze out. Although we do not 
include the Sommerfeld enhancement (SE) in our cross sec¬ 
tion computation, we show the SE for pure wino (gwino = 1) 
with the black line as a reference. Note that our wino-like (or- 
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FIG. 7: The symbols denote the scattering points in the m^-{av) parameter space that satisfy the la ranges (limits) of all the constraints 
described in the text. The bino- (red squares), wino- (orange cycles), and higgsino- (green stars) like neutralinos are defined by the neutralino 
composition fraction, g, > 0.95. In contrast, the mixed neutralino (g, < 0.95) is shown in gray diamonds. The Sommerfeld enhancement for 
pure wino (g„ino = 1) is presented here as the black line for reference. The HERD 5-yr sensitivities as shown in Fig.|5]are over-plotted. The 
left panel is for yy channel, and the right panel is for yZ channel, respectively. We would like to emphasize that in our analysis we assumed a 
non-thermal relic scenario, where the neutralino can be reproduced by late decays from other particles. However, we require that dark matter 
is not thermally overproduced in the early Universe. 


ange) points do not denote wino-wino annihilations only. In 
this region the winos are mixed with some small fraction of 
higgsinos so that the cross section can be higher than the pure 
wino annihilation cross section (black line), in particular in 
the low mass region. One can see that, even without the SE, 
HERD is able to probe the wino region with < 500 GeV 
in the yy search and < 1 TeV in the yZ search. The wino 
parameter space can be further constrained if SE is included. 
The bino region on the other hand is still hard to be probe. 
Note, however, that the points in parameter space where dark 
matter is entirely produced via thermal freeze-out are not in 
the reach of the HERD instrument. 


V. CONCLUSION 

In this work, we discuss the capability of y-ray line searches 
of HERD onboard China’s space station. Based on the de¬ 
tailed simulations of the detector performance, including the 
energy resolution, effective area and exposure, and the back¬ 
ground rejection power, the sensitivity of monochromatic 
photon detection of HERD is investigated with the maximum 
likelihood analysis. Different DM density profiles, and hence 
different optimal ROIs, are discussed. We find that the elec¬ 
trons (dominant at low energies) and Galactic diffuse y-rays 


(dominant at high energies) constitute the main contributions 
to the background. Accordingly the results are more uncer¬ 
tain at low energies (< 100 GeV) due to the different ROIs for 
the different DM density profiles. Compared with the current 
Eermi results, HERD would improve the limits by a factor of 
a few, for similar observation time. This is mainly due to the 
largely improved energy resolution of HERD (~ 1%). Com¬ 
pared to ground based lACTs such as CTA, HERD will play 
a complementary role in y-ray line searches. Eor energies be¬ 
low a few hundred GeV HERD will be more sensitive than the 
lACTs. 
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